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ABSTRACT 

The quasar in the Hubble Deep Field South (HDFS), J2233-606 (z em = 2.23) has been exhaustively 
observed by ground based telescopes and by the STIS spectrograph on board the Hubble Space Tele- 
scope (HST) at low, medium and high resolution in the spectral interval from 1120 A to 10000 A. The 
combined data give continuous coverage of the Lyman-a forest from redshift 0.9 to 2.24. This very large 
base-line represents a unique opportunity to study in detail the distribution of clouds associated with 
emitting structures in the field of the quasar and in nearby fields already observed as part of the HDFS 
campaign. Here we report the main properties obtained from the large spectroscopic dataset available 
for the Lya clouds in the intermediate redshift range 1.20 — 2.20, where our present knowledge has 
been complicated by the difficulty in producing good data. The number density is shown to be higher 
than what is expected by extrapolating the results from both lower and higher redshifts: 63 ± 8 lines 
with log Nhi > 14.0 are found (including metal systems) at < z >= 1.7, to be compared with ~ 40 
lines predicted by extrapolating from previous studies. The redshift distribution of the Lyman-a clouds 
shows a region spanning z ~ 1.383 — 1.460 (comoving size of 94 hg 5 Mpc, Q = 1) with a low density 
of absorption lines; we detect 5 lines in this region, compared with the 16 expected from an average 
density along the line of sight. The two point correlation function shows a positive signal up to scales of 
about 3 hg 5 Mpc and an amplitude that is larger for larger HI column densities. The average Doppler 
parameter is about 27 km s _1 , comparable to the mean value found at z > 3, thus casting doubts on 
the temperature evolution of the Lya clouds. 

Subject headings: cosmology: observations - quasars: absorption lines - quasars: individual: J2233-606 



1998; Outram et al. 1998; Ferguson et al. 1999) and imag- 
ing (Gardner et al. 1999). 

Here we present the interesting features shown by the 
distribution of the Lya clouds obtained by the fitting of 
the absorption lines in the redshift range 1.20 < z < 
2.20 in the high and medium resolution spectra taken 
with STIS/HST (Ferguson et al. 1999 2 ) and UCLES/AAT 
(Outram et al. 1998). The spectral resolution ranges from 
FWHM = 10 km s" 1 in the interval AA = 2670 - 3040 A, 
to 50 km s- 1 in AA = 3040 - 3530 A, and to 8.5 km s" 1 in 
AA = 3530 - 3900 A. 

2. THE DOPPLER PARAMETER AND THE COLUMN 
DENSITY REDSHIFT DISTRIBUTION 

The line parameters have been obtained using the MI- 
DAS package FITLYMAN (Fontana & Ballester 1995) in 
the spectral range 2670-3900 A through \ 2 minimization 
of Voigt profiles. At shorter wavelengths, the Lyman limit 
of the system at z ~ 1.943 absorbs a large fraction of 
the QSO flux, preventing the fit of absorption lines. The 
optical UCLES/AAT and STIS/HST spectra have been 
combined to simultaneously fit the Lyman series of the 
Lya clouds and that gives more robust results than using 

1 Based on observations made with the NASA/ESA Hubble Space Telescope by the Space Telescope Science Institute, which is operated by 

AURA, Inc., under NASA contract NAS 5-26555. 

2 See also the HDFS web site at http://www.stsci.edu/ftp/observing/hdf/hdfsouth/hdfs.html 
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1. INTRODUCTION 

A satisfactory connection between quasar absorption 
lines and galaxies close to the QSO line of sight has been 
provided so far only by relatively low redshift studies 
[z < 1.5), as galaxies at high redshifts are much harder 
to observe. For z > 2, large samples of quasar absorp- 
tion lines have been provided because at these redshifts 
the many strong UV lines are redshifted to the optical 
range, where very sensitive ground based observations can 
be obtained. However a complete picture in the whole ob- 
servable redshift range is missing because of the lack of 
information at intermediate redshifts (1.5 < z < 2). 

With the advent of the high resolution UV spectrograph 
Space Telescope Imaging Spectrograph (STIS), a detailed 
study of the intermediate redshift Lya forest and the con- 
nection to galaxy properties can now be obtained. In- 
deed, this was one of the motivations for the Hubble Deep 
Field South (HDFS). Much of the ground-based and HST 
HDFS campaign (Williams et al. 1999) has been devoted 
to the observations of the quasar J2233-606 (z em — 2.23) 
that lies in the STIS field (about 5 and 8 arcmin from the 
WFPC2 and NICMOS fields). The study of the quasar in- 
corporates both spectroscopy (Sealy et al. 1998; Savaglio 
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the Lya line alone. Only the parameters of Lya clouds 
at z < 1.60 have been obtained using the Lya absorption 
line alone. The redshift distribution of the Doppler pa- 
rameters and of the column densities are shown in Fig. 1, 
together with the instrumental Doppler width and column 
density 4<r detection limit for the Lya lines. The sam- 
ple does not include HI column densities associated with 
known metal systems, and lines close (z > 2.2) to the 
quasar redshift that are considered to be affected by the 
ionization of the quasar. The number of identified inter- 
vening metal systems is 5. Three of them have a complex 
structure (for a total of 10 Lya components) with HI col- 
umn density in each system log Nhi > 16. These clouds 
are presumably associated with galaxies and we consider 
them a different population of clouds with respect to lower 
HI column density clouds of the intergalactic medium, al- 
though many of them are presumably polluted of metals 
by nearby galaxies (Cowie & Songaila, 1998). Metal sys- 
tems have been included in the sample only to compare 
the number density redshift evolution with previous stud- 
ies at lower and higher redshifts (see section 3). The total 
number of Lya clouds excluding metal systems is then 210. 
The reported 4cr detection limit is determined by consid- 
ering the error array of the spectra and does not take into 
account the line blending effect; at z > 1.7 it is not very 
important because the line density at intermediate red- 
shifts is relatively low, but it becomes particularly strong 
at lower redshifts, due to the presence of Lyman series 
lines of the higher redshift Lya clouds that prevents de- 
tection of weak lines (Nhi < 10 13 ' 5 cm -2 ). We therefore 
consider the completeness limit in the whole redshift range 
of N HI > 10 14 cm -2 . 
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Fig. 1- The Doppler parameter (lower panel) and column 
density (upper panel) vs. redshift of the Lya forest along 
the J2233-606 line of sight. The solid and dotted lines in 
the lower panel are the instrumental Doppler widths along 
the spectra for the Lya and Ly/3 forests respectively. The 
dashed, the dotted and the solid lines in the upper panel rep- 
resent the 4a HI column density detection limit in the case 
of Doppler parameter of 40, 30 and 20 km s _1 respectively. 



Some interesting features can be noticed in the column 
density redshift distribution. There is a peak of high HI 
Lya clouds at z = 1.92 — 1.99 (~ 60 hg 5 Mpc comoving). 
A metal system with a z ~ 1.943 Lyman limit was first 
detected from the test STIS observations of the quasar. 
Therefore this peak might indicate the presence of a high 
density of galaxies at those redshifts. At lower redshifts 
(in the interval 1.383 < z < 1.460), we see a region with 
a low density of lines. The corresponding line density is 
dn/dz = 65, compared with an observed mean over the 
whole range of dn/dz = 210 lines. There are no lines with 
HI column density larger than the completeness limit of 
10 14 atoms/cm" 2 , while from the mean observed in the 
whole redshift range we would expect 4 lines. Although 
this is not statistically very significant, it is suggestive of 
the presence of a "void" of comoving size of 94 Mpc. 
The redshift measurements of galaxies from multi color 
ground observations and from the very deep STIS images 
of the field will probably confirm if this is real. The lack 
of Lya lines can also be caused by a strong ionization of 
the clouds by the local UV radiation field, for instance in 
a quasar environment. 

It is worth noting the decrement of lines going from the 
QSO redshift to ~ 1.38, shown by the histogram of the 
entire sample (Fig. 2). This effect disappears when select- 
ing strong lines, and so it may be due to the detection 
limit variation along the spectrum. For smaller redshifts, 
the number of lines increases again. We notice that at 
z ~ 1.335, a CIV system has been tentatively identified 
(Ferguson et al. 1999) and a quasar at a distance of about 
44.5" from the line of sight (~ 300 h^ 5 kpc) has been found 
from the ground by EMMI/NTT observations (Tresse et 
al. 1998). This might also be an indication of an overden- 
sity of objects at that redshift. 
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Fig. 2- Histogram of the HI column density of Lya clouds 
as a function of redshift for different thresholds. Error bars 
are the square-root of the number in each bin. 
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3. THE NUMBER DENSITY EVOLUTION AT Z = 1.20 - 2.20 

The number density evolution of the Lya forest, de- 
scribed by a power law of the type dn/dz oc (1 + z) 1 , has 
been studied at high rcdshift (1.7 < z < 4.1) using sam- 
ples of high resolution data by Giallongo et al. (1996) and 
Kim et al. (1997). For lines with \ogN HI > 14.0, dn/dz 
shows a fast evolution, with 7 ~ 3.6. The analysis at 
low redshifts (Weymann et al. 1998), based on low resolu- 
tion spectroscopy of the Quasar Absorption Line Survey 
(QALS), gives in the range 0.0 < z < 1.5 a much flat- 
ter dn/dz, consistent with weak negative evolution (for 
Q a = 1) 7 = 0.16 ± 0.16. We note that in the two redshift 
regimes a variation of 7 with column density threshold has 
been found (larger for increasing threshold) . This would go 
the opposite direction than what is shown in Fig. 2 where 
7 is smaller for higher threshold. The discrepancy can be 
explained if a large number of lines with log Nhi < 14 is 
present, but not detected at z = 1.4 — 1.8. 
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Fig. 3- Number density evolution of the Lya clouds with 
log Nhi > 14.0 from 2 = to z = 4. Filled symbols are for 
samples that include metal systems, open symbols do not. 
In the sample of Kim et al., only lines with log Nhi < 16.0 
are included. The low redshift line is the fit with 7 = 0.16 
obtained by a sample of Lya clouds with equivalent width 
EW > 0.24 A (Weymann et al. 1998) which corresponds to 
log Nhi > 14.0 for b ~ 26 km s -1 . The high redshift line is 
the fit with 7 = 3.62. 

The problem in comparing high and low redshift data 
stems from the different techniques used to measure col- 
umn densities. At high redshifts, high resolution obser- 
vations allow line fitting of Voigt profiles, and those give 
more robust results than the curve-of-growth technique 
used to estimate HI from equivalent width measurements 
of low resolution data at low redshifts. On the other hand, 
the strong evolution of line density makes the line deblend- 
ing of complex structures rather difficult at high redshifts. 
The Lya forest of J2233-606 lies between the two red- 
shift regimes and so it is extremely useful to evaluate the 
correct connection between the two. The results for the 
two samples of Lya clouds with and without metal lines 
are shown in Fig. 3 for lines with log Nhi > 14.0. The 



data from J2233-606 have been binned into three redshift 
ranges Az = 1.201 - 1.535, 1.535 - 1.870 and 1.870 - 2.204 
corresponding roughly to the STIS high resolution and 
medium resolution intervals and to the UCLES high res- 
olution interval. The last two bins give a number density 
of lines that would better match the high and low redshift 
observations if metal systems are excluded. The first bin 
gives a larger number of lines, 66 ± 14 including metal sys- 
tems, compared with 40 ± 4 lines given the QALS best 
fit (Weymann et al. 1998). That deviation is not particu- 
larly significant, and could indicate a peculiarity along the 
line of sight to J2233-606. However, a similar excess at 
1.505 < z < 1.688 (63 ± 14) is also shown by the quasar 
UM18, which is part of the FOS sample of the QALS, but 
not included in the fitting of the number density evolution. 
The excess found in J2233-606 can be due to the presence 
of a cluster of lines at z ~ 1.335, as suggested in the pre- 
vious section. Other possible explanations would be the 
blending effect by which nearby lines with column den- 
sity below the threshold would artificially appear in the 
limited signal-to-noise spectrum as stronger single lines 
with HI above the threshold, or the presence of a small 
number (~ 5 would be enough) of unidentified metal lines 
interpreted as Lya lines. This latter possibility looks very 
unlikely because even in clusters, metal lines typically ap- 
pear very narrow at high resolution. 
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Fig. 4- Two point correlation function of the Lya forest in 
J2233-606 for lines with log Nhi > 13.0 and 13.8. The Pois- 
son error is also given as dashed and solid lines respectively. 

4. THE CLUSTERING PROPERTIES 

A positive signal in the clustering properties of the 
Lya forest has been found thanks to the advent of high res- 
olution spectroscopy. The investigation has been mainly 
based on the analysis of the two point correlation func- 
tion (Cristiani et al. 1997; Kim et al. 1997) but also on 
that of the power spectrum of mass density fluctuations 
(Hui et al. 1997; Croft 1998; Amendola & Savaglio 1998). 
Results have shown that the clustering is present up to 
Aw ~ 300 km s _1 and that the amplitude increases with 
increasing column density threshold of the lines. Most 
interesting for models of large scale structure formation 
is the redshift evolution found, showing larger values at 
smaller redshifts. Although the analysis of the two point 
correlation function is uncertain when the sample is too 
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small, we report the result on the forest of J2233-606 in 
Fig. 4 for two column density thresholds. The signal is sig- 
nificant in the first bin (D ~ 0.8 Mpc) and it is higher 
at the higher HI threshold at about the 6 and 4 a level for 
logiVffj > f3.0 and f3.8 respectively. In a previous study 
by Cristiani et al. (1997), the two point correlation func- 
tion for logNni > 13.8 at Av = 100 km s _1 has shown a 
rcdshift evolution, being £i o — 0.20 at < z >= 3.85 (con- 
sistent with no clustering), £i o — 0.75 at < z >= 3.40, 
and ^ioo ^ 0.85 at < z >= 2.40. The line of sight of 
J2233-606 shows £ioo — 1-3 at < z >= 1.7 and confirms 
this evolution. 

5. THE DOPPLER PARAMETER DISTRIBUTION 

Thanks to the high resolution observations of J2233- 
606, the analysis of the distribution of the Doppler pa- 
rameter of the Lya forest can be performed for the first 
time at these low redshifts. At higher rcdshifts Kim et al. 
(1997) have found that the median value of the Doppler 
parameter is b ~ 30 km s _1 at a rcdshift of z ~ 3.7, and 
appeared to increase to b ~ 35 — 40 km s _1 at a rcdshift 
of z ~ 2.3 for lines with 13.8 < log Nhi < 16, which has 
been interpreted as an increase of the temperature or kine- 
matic broadening of the clouds. The Doppler distribution 
in J2233-606 is shown in Fig. 5 for the whole sample and 
for the subsample of lines (60% of the total) with smaller 
uncertainties (cr(b) < 8 km s _1 and a (log Nhi) < 0.5). 
The fit of the two Gaussians (in the interval < b < 60 
km s _1 ) gives in the two cases a mean value of 27 and 
26 km s -1 respectively (x 2 is 2.4 and 1.5), and these are 
slightly different than the median value (29 and 25 km s _1 ) 
in the whole b range. The median values are larger if lines 
with 13.8 < log Nhi < 16 are selected, being 31 km s _1 in 
the two cases. Thus, the median value we find over the 
rcdshift range 1.20 < z < 2.20 is considerably below the 
value of b ~ 40 km s _1 expected on the basis of the redshift 
evolution suggested by Kim et al. (1997). 

6. CONCLUSIONS 

In contrast to the empty field used for the HDF North 
program, the HDFS has been selected close to a high 
redshift quasar, J2233-606 (z em = 2.23). The QSO is 
centered in the STIS field, located approximately 5 and 
8 arcmin from the WFPC2 and NICMOS fields respec- 
tively. The Lya forest of J2233-606 has been observed 
at high and medium resolution (FWHM = 10, 50 and 8.5 
km s~ 4 in the wavelength intervals AA = 2670 — 3040 
A, AA = 3040 - 3530 A and AA = 3530 - 3900 A) us- 
ing STIS/HST (Ferguson et al. 1999) and UCLES/AAT 
(Outram et al. 1998), giving the opportunity to analyze 
the properties of the intergalactic medium in the large 
rcdshift interval z — 1.20 — 2.20. This is the first time 
that the Lya forest is studied in detail for z 2. More 
than 200 Lya cloud redshifts have been found and 63 have 
HI column density in excess of 10 14 atoms/cm -2 , 11 of 
which have associated metal absorption. The redshift dis- 
tribution of the lines has shown a deficiency of lines (5 
are found over a mean expected of 16) in the interval 
1.383 < z < 1.460, corresponding to a comoving size of 94 
Mpc (fi = 1). The region around the metal system 



at z = 1.942, also responsible of the Lyman limit observed 
for A < 2700 A, shows an excess of strong lines in the inter- 
val 1.92 < z < 1.99, indicating a clustering of Lya clouds. 
An excess of lines is also found for 1.20 < z < 1.38. The 
evolution of the number density in the total sample shows 
a rapid decrease of lines with redshift. This is partially 
due to the blending effect of low HI column density lines, 
which is stronger at lower redshifts. Indeed it tends to 
disappear when selecting only strong lines (Nhi > 10 14 
cm -2 ). On average the number of strong Lya clouds per 
unit redshift is higher than the extrapolation from studies 
at lower and higher rcdshifts: 63 ± 8 at < z >= 1.702 
including metal systems are found to be compared with 
41 and 27 expected from the low and the high rcdshift 
predictions. A larger sample of lines will confirm if this 
line density is typical at that redshift or if it is peculiar 
to the J2233-606 line of sight. The two point correlation 
function has shown a positive signal up to a scale of the 
order of 3 Mpc. The amplitude found seems to con- 
firm the evolution with redshift, being larger at lower z. 
The mean Doppler parameter of the Lya forest is around 
26 — 27 km s _1 , corresponding to a temperature, in the 
case of thermal broadening, of ~ 4.3 x 10 4 K. The me- 
dian Doppler parameter evolution reported by Kim et al. 
(1997) is not confirmed by this Lya forest. 

It is a pleasure to thank L. Amcndola, B. Carswell, S. 
Casertano, G. Ganis, B. Jannuzi, M. Livio, P. Outram and 
R. Weymann for useful suggestions. 
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Fig. 5- Doppler parameter histograms of Lya clouds in 
J2233-606 in the whole sample (upper panel) or in the sub- 
sample of lines with la error in the column density and in 
b smaller than 0.5 dex and 8 km s _1 respectively (124 lines 
over a total of 210). The fit of a Gaussian for b < 60 km s _1 is 
shown as well. Error bars are the square-root of the number 
of lines in each bin. 



5 



REFERENCES 



Amendola L., Savaglio S., 1998, MNRAS, submitted 

Cristiani S., D'Odorico, S., D'Odorico, V., Fontana, A., Giallongo 

E., Savaglio S., 1997, MNRAS, 285, 209 
Cowie L. L., Songaila A., 1998, Nat, 394, 44 

Croft R. A. C, Weinberg Pi H Ppftini M.., Hrtnquist L., Katz N., 

1998, ApJ, submitted, astro-ph/9809401 
Ferguson H. C. et al., lyyy, in preparation 
Gardner J. P., 1999, in preparation 

Giallongo E., Cristiani, S., D'Odorico S., Fontana A., Savaglio S., 

1996, ApJ, 466, 46 
Fontana A., Ballester P., 1995, The ESQ Messenger, 80, 37 



astro-ph/9809404 



Hui L. Gnedin N. Y., Zhang Y., 1997, ApJ, 486, 599 
Kim T-S, Hu E. M., Cowie L. L., Songaila A., 1997, AJ, 114, 1 
Outram P. J., Boyle B. J., Carswell R. F., Hewqtt P, C, William 

R. E., Norris R. P., 1998, MNRAS, submitted 
Savaglio S., 1998, AJ, 116, 1055 
Sealey K. M., Drinkwater M. J., Webb J. K., 1998, ApJL, 499, L135 
TVoqqc T. ncmnofolrl M., Petitjean P. et al., 1998, A&A, submitted, 

astro-ph/981224C 
Wcymann R. J., Jannuzi B. T., Lu L., et al., 1998, ApJ., 506, 1 
Williams R. E. et al., 1999, in preparation 



